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We study link-diluted ± J Ising spin-glass models on the hierarchical lattice and on a three-dimensional lattice 
close to the percolation threshold. We show that previously computed zero temperature fixed points are unstable 
with respect to temperature perturbations and do not belong to any critical line in the dilution-temperature 
plane. We discuss implications of the presence of such spurious unstable fixed points on the use of optimization 
algorithms, and we show how entropic effects should be taken into account to obtain the right physical behavior 
and critical points. 
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Introduction. — Frustrated systems may have very complex 
free-energy landscapes at low temperature T which in turn 
give rise to very peculiar thermodynamic al properties 
IslH]- It was recently shown in Ref. ijH] for the case of 
the two-dimensional Edwards-Anderson (EA) model that in 
such situations T = computations as the ones of Ref. |@] 
may produce misleading results. In Ref. [7| it was shown 
that there are ways to improve on the results of Ref. [6] us- 
ing T — methods; however the r/ exponent is still found to 
be non-universal; this could be due to entropic effects, which 
are neglected in Ref. iQ] ■ In this Letter we study entropic ef- 
fects in the very low temperature regime of disordered mod- 
els. We consider diluted spin glasses (SG) with discrete cou- 
pling distributions at very low temperature and show that pre- 
vious computations of critical points were incorrect JES- We 
explain why and we show how to compute the right critical 
points: considering first-order corrections in T allows us to 
treat the entropic contribution to the free energy correctly in 
the limit when T goes to zero. A similar idea has been applied 
also to models defined on random graphs lllOll . Here we find 
that, on hierarchical as well as three-dimensional (3D) lattices, 
the SG phase persists in a region of higher dilution than calcu- 
lations done exactly at T = would predict. In other words, 
we find SG ordering induced solely by entropic effects which 
is a phenomenon somewhat reminiscent of Villain's "order by 
disorder" (|Ii|ljl]. 

We consider two SG models: one defined on the hierarchi- 



cal lattice III2II . where the Migdal-Kadanoff approximation is 
correct, and the other one on a 3D lattice. In the former case, 
we are able to compute the distribution of effective couplings 
and of the free energies allowing us to determine exactly the 
location of the critical points and the whole phase diagram. 
For the 3D case, we have to resort to numerical simulations: 
using state-of-the-art Monte Carlo (MC) methods for disor- 
dered systems, we provide evidence that previously computed 
critical lines are likely to be incorrect, and by a numerical 
study of percolation properties we give bounds on the loca- 
tion of the right critical points. 
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FIG. 1: Phase diagram of the diluted ±J spin glass in the tempera- 
ture versus bond density plane, where pc denotes the geometric per- 
colation threshold, psG the percolation threshold for spin-glass order 
and p* the percolation threshold for spin-glass order determined at 
T — 0. Inset: Elementary step for the construction of a hierarchical 
lattice with 6 = 3 and s = 2. 



Model. — We consider an Ising SG model defined by the 
Hamiltonian 



(1) 



with Ising spins ai = ±1 and the sum running over all the 
nearest neighbors pairs on the lattice. The couplings Jy 
are quenched, independent and identically distributed random 
variables extracted from the distribution 11311 



Po(./) = (1 - p)SiJ) + ^ 5{J - 1) + 5{J + 1) 



(2) 



We are interested in the T — Q critical point separating the 
paramagnetic phase {p < pso) from the SG phase (p > psg)- 
We consider two lattices: the hierarchical lattice and the 3D 
simple cubic lattice for which the Hamiltonian in Eq. ([T]i cor- 
responds to the EA model 1 14]. 

Hierarchical lattice. — The hierarchical lattice of G gener- 
ations is obtained by applying the construction shown in the 



inset of Fig. [T| to all the links of the G ~ 1 generations lat- 
tice (the G = lattice being a single link connecting two 
vertices). Such a construction is defined by two parameters: 
the number b of parallel branches, made of s bonds in se- 
ries each. The effective dimension of the model is defined by 
d = 1 + ln(6) / ln(s) and it is known that the lower critical di- 
mension for the SG transition is close to di ~ 2.5 ll6i [Tsi [l6ll . 
Working with d > diwe have Tsg > for p — 1 and a criti- 
cal line in the {p, T) plane (see Fig.[T]). We present results for 
6 = 3 and s = 2 (i.e., d = 2.585...), but we have checked 
their validity also for other choices of the parameters. 

A model defined on the hierarchical lattice can be solved 
exactly by recursive decimation: starting from a lattice of G 
generations and summing over the spins introduced in the last 
generation one can easily obtain a lattice of G — 1 genera- 
tions with renormalized couplings. At each step of this renor- 
malization procedure the couplings undergo two elementary 
transformations: (i) each s-tuple of couplings in series pro- 
duces an effective coupling J and (ii) each 6-tuple of effec- 
tive couplings in parallel linking the same pair of variables is 
summed together, giving the renormalized coupling. In the 
T ^ limit, one renormalization step can be written as 

J(") 4 sign(j(") ..... J(")) mind . . . , I ^i"^!), (3) 

j^n+l) d Jin) ^ ^ ^ ^ ^ jin)^ (4) 

where equalities hold in distribution sense and the index (n) 
stays for the number of renormalization steps. Applying these 
functional recursion equations, the initial distribution Pq of 
couplings flows to the fixed-point one Poo- The critical den- 
sity of links p* is defined such that P^ is nontrivial, i.e., 
different from the paramagnetic one, 6{J), and the SG one, 
[S{J - oo) + S{J + oo)]/2. 

In Ref. 101 Bray and Feng computed the distribution P^ 
under the assumption that the renormalized couplings remain 
integer, which is indeed a possible solution to Eqs. Q and 
Actually in Ref. iQ] only three values for the couplings (0, ±1) 
were considered; the extension to a symmetric distribution in- 
volving all the integer values between —AI and AI is straight- 
forward, and leads in the limit of large M to p* ~ 0.465. In 
the inset of Fig.|2]we show the variance of J(") as a function 
of n for various values of p. This critical point p* has been 
considered up to now the boundary between the paramagnetic 
and the SG phases. Nevertheless, we find that in the plane 
{p, T) the point (p*, 0) does not belong to any critical line; 
i.e., it is an isolated point (see Fig.[T]l, and for this reason it 
is irrelevant for the thermodynamics at any positive T. Even 
at T = it does not correspond to the border between two 
different phases, i.e., p* ^ psg' 

In order to prove that the fixed point [p* , 0) is isolated as in 
Fig. [U we study coupling renormalization at infinitesimal T. 
We start from the positive T rule for decimating s = 2 bonds 
in series (for s > 2 the same transformation must be applied 
more than once) 

tanh(/3J) = tanh(/3Ji) tanh(/3J2), (5) 
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FIG. 2: The probability that the variance of the renormalized "en- 
tropic" couplings K goes to zero as a function of bond density p, for 
different sizes A'^ of the population evolved numerically. The inset 
shows the behavior of the variance of the renormalized "energetic" 
couplings I, obtained analytically in the N —> oo limit. 

with (3 = 1/T. Similarly to iH, for very low T, we 
rewrite J = TK in case J is vanishing for T ^ 0, and 
J = sign(/)(|/| - TK) if limT^o J = / 7^ 0. The choice 
in the sign is dictated by the fact that thermal fluctuations de- 
crease the coupling intensity. We have checked that higher or- 
der corrections in temperature are unnecessary in the T ^ 
limit. Plugging these new variables in Eq. (|5]l we find that 
in the T ^ limit the "energetic" component / gets renor- 
malized following Eq. Q, while the "entropic" component K 
follows 

K^Ki if|/i|<|/2| (6) 

tanh(7^) = tanh(Ki) tanh(if2) if h = h = (7) 
exp(2K) = exp(2ii'i) + exp(2is:2) if h = h 7^ 0. (8) 

Before renormalization, the original couplings are T-indepen- 
dent and distributed according to Pq; thus we start with / e 
{0, ±1} and K — 0. The recursive equation for / is the one 
already studied above and gives a coupling flowing to for 
p < p* (see inset of Fig. |2j. Nevertheless, even if / = for 
all the renormalized couplings, the entropic correction may 
be non-null, thus giving a nonzero correlation of the order of 
tanh(/3 J) = tanh(i^) ^ 0, even in the T ^ Hmit. Please 
note that if the temperature was set to zero at the beginning of 
the calculation such a nonzero correlation would not be found. 

Under renormalization the variance of entropic couplings 
goes to zero or diverges, depending on the value of p. In the 
main panel of Fig. |2] we show the probability that {K'^)i=q 
flows to zero, where N is the size of the population used 
for simulating the decimation procedure. It clearly shows 
the existence of a critical value psc — 0.4588 such that for 
P < PSG correlations decay to zero at large distances, while 
for p > psG there are long-range correlations (of SG type, 
since (A'} ~ 0). 

On hierarchical lattices, the inequality psG < P* holds in 
general, and implies that for T = and p G {psG , P* ) the long 
range-order is induced solely by entropic effects. That is the 
two ground states obtained by fixing the outermost spins of the 
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hierarchical lattice in a parallel or anti-parallel way have ex- 
actly the same energy, but very different entropies; long-range 
correlations from the ground state with the largest entropy will 
dominate ensemble averaged correlations. Such long-range 
correlations arise because when summing two parallel effec- 
tive bonds such that Ii + 12 = it may be that Ki + K2 7^ 0, 
i.e., there is no perfect cancellation and an effective coupling 
of intensity 0{T) persists. In the T ^ limit these couplings 
induce nonzero correlations, which may eventually increase 
under renormalization (as for p > psc)- It is worth noticing 
that this mechanism for generating 0{T) couplings may work 
perfectly well also when applying a decimation procedure on 
regular lattices [9]. 

A similar argument also tells us that on the very first steps of 
decimation, when K takes few and discrete values (for large 
n, K becomes dense on R), exact cancellations, both in / 
and K, take place. If the original model had a bond density 
slightly larger than the percolating density pc (for 6 = 3 and 
s — 2, pc ~ 0.389391...) these exact cancellations would 
produce a renormalized lattice which is below the percolating 
point, preventing any long-range correlation from arising. So 
also the inequality pc < psG must hold in general for any 
frustrated model with discrete couplings on the hierarchical 
lattice. 

Regarding the issue of the universality, we checked that 
starting from any point on the critical line extending from 
(psG,0) to (IjTsg), see Fig. [T] the same fixed-point cou- 
pling distribution is obtained; note that, at a positive tempera- 
ture, couplings J/T get renormalized exactly by Eq. (|7]l. This 
fixed-point distribution has nothing to do with the integer- 
valued one obtained from {p* , 0). 

3D simple cubic lattice. — The phase diagram shown in 
Fig. [T] is exact for the SG model defined by Eq. ([T]) with 
coupling distribution from Eq. (|2|i on the hierarchical lattice 
(or equivalently in the Migdal-Kadanoff approximation). We 
present now evidence that the main qualitative features of that 
phase diagram are preserved when a 3D lattice is considered. 
With this in mind, we show by MC simulations and by perco- 
lation arguments that (i) a SG transition takes place at finite 
temperature for p < p*, thus implying pso < P* and that (ii) 
there exists a percolating phase with no long-range order, thus 
implying Pc < PsG with pc = 0.2488126(5) [18]. 

The MC technique we use to simulate the 3D link-diluted 
± J EA model close to the percolation threshold is a combi- 
nation of the replica cluster update moves 1 19 1 embedded in 
parallel tempering |20] as described in Ref. 112 III in alternation 
with standard Swendsen-Wang f23\ cluster updates on each 
of the replicas. The additional Swendsen-Wang updates are 
very efficient close to percolation and speed up the simulation 
considerably |23|]. This is an important issue because close 
to percolation we encounter noticeable finite-size effects, and 
therefore relatively large lattices have to be simulated to see 
a good signal for a SG transition below p*. We use the value 
p* — 0.272(1) determined by Boettcher 191] from the behavior 
of the domain-wall defect energies of the ground states using 
an optimization algorithm ll24ll . 
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FIG. 3: (Color online). £^{L,T)/L as a function of temperature for 
several system sizes L and different bond densities p. 



To check for the presence of a SG transition we measure 
the second-moment correlation length ^(L, T) on a lattice of 
size L defined as 112511 
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(9) 



2sin(|kmin|/2) 
where the wave-vector-dependent SG susceptibility is 



XSG( 



and k,„in = (0, 0, 2tt/L) is the smallest nonzero wave vector 
allowed by periodic boundary conditions. We denote thermal 
averages by (•) and disorder averages by ~. 

In Fig.[3]we display the behavior of ^{L, T)/L for four dif- 
ferent values of p. In the upper left-hand panel we show the 
result below the percolation threshold {p = 0.24 < pc) for 
which there is no SG transition as can be seen from the ab- 
sence of a crossing of £^{L,T)/L for different lattice sizes. 
The upper right-hand panel shows the situation for pc < p = 
0.2625 < p* where again the curves do not cross and therefore 
there is no sign for a SG transition. For p — 0.268, slightly 
below p*, we find — amongst noticeable finite-size effects 
— that the ^(L, T)/L for different lattice sizes do cross in- 
dicating the presence of a SG transition. The crossing point 
is moving slightly towards a smaller effective Tc as L is in- 
creased, but this is also the case for p = 0.2725 ^ p* and 
larger values of p. The fact that the crossing happens at a value 
of ^(L, Tc) /L « 0.6, with a tendency to increase for larger L, 
is an important evidence for a SG transition since this value is 
a renormalization group invariant quantity and was found to 
lie between 0.60 and 0.65 in different recent studies of the 3D 
SG transition lllllliei |27ll . 
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the crossover length even larger This argument suggests that 
any numerical method at zero temperature is plagued by the 
unphysical fixed point up to unreachable length scales! 
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IST-034952 (GENNETEC). 
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FIG. 4: (Color online). Percolation probability in a 3D ± J diluted 
spin-glass model, reduced with T = 0^ rules. 

To show that indeed pc < psG holds in 3D, we apply to 
the Hamiltonian in Eq. ([TJ the reduction rules for weakly con- 
nected lattice sites (i.e., with degree 1 or 2) as described in 
Ref. i^. Instead of applying these rules exactly at T = we 
again keep terms of order T for the resulting effective cou- 
plings leading to a more connected remaining graph than the 
same procedure at T = would produce. We apply the reduc- 
tion rules recursively until the graph cannot be reduced any 
further. At this point we determine whether there exists at 
least one path, along which the remaining effective couplings 
are nonzero, that connects the z = and z = L plane of the 3D 
lattice with periodic boundary conditions in the x and y and 
open boundary conditions in the z direction. The presence of 
such a path is a necessary, but not sufficient condition for a 
SG phase. In Fig.|4]we show the probability that such a path 
exists as a function of p for different sizes L. It clearly shows 
that exact cancellations also in 3D lead to pc <psG- 

Discussion. — An important comment is in order: opti- 
mization methods that compute ground state energies as in 
Ref. 1 9] are typically sensible to the "unphysical" isolated 
critical point at p* rather than to the "physical" critical point 
at psG- To detect the right critical point these optimization 
methods should be modified to include entropic effects. In 
Ref. |28] an attempt in that direction was made, leading to a 
change in physical observables. The decimation method used 
in Ref. [9] can be easily extended to consider first-order cor- 
rections in temperature, thus leading to a more connected lat- 
tice after decimation. 

However, even if a correct procedure is employed, taking 
into account entropic effects, the presence of the spurious 
fixed point at p* may lead to very strong crossover effects. 
For example for p £ {psG ,P*), as long as energetic couplings 
are present, they dominate the behavior of any observable; the 
right physical behavior, given by entropic couplings, can be 
measured only at length scales where energetic couplings have 
disappeared, and these scales may be extremely large. In the 
b — 3 and s = 2 hierarchical lattice such a length scale is 
always larger than s^^ = 32768, roughly. Increasing the di- 
mensionality we find that psG and p* get closer, thus making 
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